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ABSTRACT 


Fluctuating  temperature  structures  in  the  atmosphere 
induce  phase  perturbations  in  a  propagating  laser  beam.  These 
turbulent  conditions  occur  throughout  the  atmosphere  and  cause 
the  laser  beam  to  spread  and  alter  its  centroid.  There  are 
several  methods  to  measure  the  parameters  of  optical 
turbulence  in  the  atmosphere,  but  few  that  will  determine  them 
as  a  function  of  altitude  at  all  levels.  One  method  of 
measuring  altitude  profiles  of  turbulence  is  with  a 
temperature  probe  launched  via  a  balloon  system. 

This  thesis  involves  the  development  of  a  differential 
temperature  probe  sensor  to  measure  the  temperature 
fluctuations  at  all  altitudes  in  the  atmosphere.  In  addition, 
it  investigates  the  effect  of  solar  heating  on  the  probes  in 
the  atmosphere  and  the  subsequent  effects  on  the  measurements. 
A  validation  of  the  probe  system  was  made  by  a  comparison  test 
with  an  acoustic  echosounder  developed  earlier.  In  addition 
to  validating  the  probe  system,  the  absolute  analysis  of 
the  echosounder  was  confirmed. 
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I .  INTRODUCTION 


Atmospheric  conditions  will  cause  severe  degradation 
along  the  optical  path  of  a  ground  to  space  weapon  or 
surveillance  laser.  There  are  several  causes  for  this 
degradation,  they  are  1)  absorption  and  scattering  by  aerosols 
such  as  rain  and  clouds,  2)  distortion  by  thermal  blooming 
and  3)  distortion  by  atmospheric  turbulence  [Ref.  1] . 
Absorption  and  scattering  can  be  controlled  by  varying  the 
wavelengths  of  the  laser  and  having  multiple  sites  to  insure 
at  least  one  has  a  cloud  free  line  of  sight.  Thermal  blooming 
is  the  heating  of  the  medium,  through  which  the  laser  beam 
propagates.  It  is  due  to  absorption  of  the  radiation  by 
molecules  and  aerosols  and  the  consequent  distortion  of  the 
beam  due  to  density  reductions  brought  on  by  the  heating. 
Choosing  approximate  wavelengths  of  the  laser  which  have  low 
atmospheric  absorption  in  the  atmosphere  reduces  this  effect. 
Atmospheric  turbulence  is  difficult  to  deal  with  since  there 
is  no  way  to  avoid  it. 

The  major  effect  of  turbulence  on  an  optical  beam  is  the 
limitation  of  the  mutual  coherence  lengths.  For  example,  an 
average  ground  to  space  coherence  length  for  the  atmosphere 
is  on  the  order  of  5  cm  therefore  a  ground  based  laser  having 
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a  4  meter  diameter  mirror  will  deliver  less  then  1/1000  of  its 
original  power  onto  a  target.  Adaptive  optics  provides  a  means 
for  reducing  these  turbulent  effects.  It  corrects  for  the 
effects  of  turbulence  by  altering  the  wavefront 
characteristics  of  a  beam  [Ref.  2]  using  deformable  mirrors 
or  nonlinear  optical  materials  such  as  Barium  Titanate. 

Measurements  of  the  turbulence  from  the  surface  to  an 
altitude  where  the  turbulence  is  no  longer  significant 
(approximately  30  km.)  are  needed  for  several  reasons.  The 
most  important  is  to  be  able  to  characterize  the  turbulent 
profiles  of  the  atmosphere  at  different  locations  in  order  to 
select  the  best  site  for  a  ground  based  optical  or 
surveillance  system.  Each  site  will  have  a  different 
turbulence  profile  since  it  depends  on  the  local  geography  as 
well  as  the  upper  atmospheric  conditions  controlled  by  general 
meteorological  patterns.  Another  reason  for  the  measurements 
is  to  determine  vertical  distribution  of  the  turbulence 
parameters,  such  as  the  coherence  length,  which  affect  the 
design  of  adaptive  optics  systems,  or  the  signal  processing 
transformations  in  a  surveillance  system. 

There  are  several  instruments  for  measuring  turbulence 
parameters,  some  of  which  measure  the  index  of  refraction 
structure  parameter  and  others  that  measure  the  temperature 
structure  parameter.  Some  of  the  methods  include  analysis  of 
star  trails  on  photographic  emulsions  [Ref.  3],  an 


2 


isoplanometer  which  measures  the  isoplanatic  angle  through 
stellar  scintillation  [Ref.  4],  and  a  Modulation  Transfer 
Function  (MTF)  device  for  determining  the  coherence  length 
[Ref.  5].  These  devices  measure  important  integrated 
parameters  but  they  cannot  measure  the  vertical  profile  of 
turbulence.  An  acoustic  echosounder  which  is  similar  in  design 
and  construction  to  a  SONAR  [Ref.  6],  measures  a  vertical 
profile  but  is  usually  limited  to  several  hundred  meters 
range.  Greater  vertical  resolution  occurs  at  the  expense  of 
decreased  maximum  range. 

In  order  to  get  a  measurement  of  the  vertical  profile  of 
turbulence  to  20-50  km  a  device  must  be  able  to  be  launched 
on  a  sounding  balloon  or  an  aircraft.  An  example  of  this  type 
of  device  is  the  thermosonde  originally  designed  by  GTE 
Sylvania  and  revised  by  the  Air  Force  Weapons  Laboratory  and 
Tri-Con  Associates  Inc.  and  built  and  used  by  the  Air  Force 
Geophysics  Laboratory  [Ref,  7]. 

This  thesis  is  an  attempt  to  design  and  build  a 
temperature  sensing  probe  system  to  measure  the  vertical 
profile  of  the  temperature  structure  parameter  which  is  a 
measure  of  the  vertical  profile  of  turbulence.  Although  such 
systems  exist,  they  are  1)  expensive,  greater  than  $2000  per 
launch  and  2)  require  extensive  calibration.  The  purpose  of 
this  thesis  was  to  develop  an  inexpensive  device  (the  system 
cost  is  approximately  $150  per  launch)  which  is  simple  to 
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operate  (it  is  a  self-calibrating  device).  Additionally  the 
effects  of  solar  heating  of  the  probes  are  investigated. 

The  results  of  a  comparison  with  a  well  developed 
acoustic  echosounder  indicates  the  system  will  be  effective 
in  measuring  turbulence  in  the  atmosphere.  The  studies  of  the 
solar  heating  effects  indicate  the  only  source  of  error  due 
to  solar  heating  would  be  from  the  sun/shade  effect  on  the  two 
probes  and  if  that  is  corrected  for,  it  will  be  accurate  up 
to  30  km. 
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II.  BACKaROUMD 


A.  TURBULENCE  PARAKETERS 

Small  temperature  variations  carried  by  the  turbulent 
velocity  field  in  the  atmosphere  produce  small  phase 
perturbations  in  an  optical  plane  wave  propagating  through  it. 
These  perturbations  randomly  distort  and  convolve  the  phase 
of  a  plane  wave.  There  are  three  atmospheric  parameters  which 
must  be  determined  prior  to  any  attempt  made  at  compensating 
for  these  atmospheric  distortions.  These  parameters  are  the 
refractive  turbulence  structure  parameter,  ,  the  spatial 
coherence  length  of  the  atmosphere,  ,  and  the  isoplanatic 
angle  9^. 

The  most  important  of  these  parameters  is  Tatarski 
[Ref.  3]  states  that  one  way  to  deal  with  a  non-stationary 
problem,  which  includes  all  atmospheric  parameters,  is  to 
define  a  function  in  terms  of  a  difference,  he  than  defines 
the  structure  function  for  index  of  refraction  as, 

=<1N(?2)-N(rj)l2>  _  ,,, 

where  <  >  denotes  an  ensemble  average.  If  we  assume  the 
atmosphere  to  be  homogeneous  and  isotropic  over  small  regions 
the  structure  function  can  be  rewritten  as. 
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(2) 


D,(r)  =  <[N(r2)-N(ri)  ]S, 

where  r  is  r,-r2.  By  dimensional  analysis  Kolomogorov  [Ref.  3] 
showed  that  the  structure  function  has  an  dependency. 
Consequently  is  proportional  to  a  constant  times  r^''^, 

Dn  =  (3) 

is  the  refractive  turbulence  structure  parameter,  a  mean- 
square  statistical  average  of  the  difference  in  the  index  of 
refraction  between  two  points  separated  by  r^2» 

=  <(N2-Ni)  (4) 

The  normalization  extends  from  an  inner  scale  1^^,  on  the 
order  of  millimeters,  to  an  outer  scale  of  L^,  on  the  order  of 
meters  [Ref.  3].  These  fluctuations  in  the  index  of  refraction 
arise  from  variations  in  density  caused  by  temperature 
fluctuations  in  the  turbulent  velocity  field. 

These  density  variations  in  the  atmosphere  alter  the 
phase  of  an  optical  beam  being  propagated  through  it.  The 
Optical  Transfer  Function  (OTF)  characterizes  the  integrated 
phase  perturbations  of  an  optical  beam.  It  is  a  measure  of  the 
correlation  of  the  electric  fields  of  the  optical  beam 
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perpendicular  to  the  direction  of  propagation.  Although  the 
atmosphere  is  not  homogeneous  or  isotropic,  Tatarski  [Ref.  8] 
postulates  the  idea  of  local  homogeneity  and  isotropy,  which 
states  that  over  some  region  R,  comparable  to  the  outer  scale 
length  L^,  we  can  assume  the  atmospheric  random  variables  are 
homogeneous  and  isotropic.  The  modulus  of  the  OTF  is  the 
atmospheric  Modulation  Transfer  Function  (MTF) .  Fried  [Ref. 
9]  introduces  the  parameter  r^  to  characterize  the  MTF.  It  is 
related  to  the  refractive  turbulence  structure  parameter  by. 


where  r^  is  the  spatial  coherence  length,  k  is  the  wave 
number  {Ir/x),  and  is  the  refractive  turbulence  structure 
parameter  along  the  optical  path  of  length  L  [Ref.  5]. 

The  other  measure  of  spatial  coherence  in  the  atmosphere 
is  the  isoplanatic  angle  9^,  It  is  similar  to  r^  in  that  it 
is  the  dependence  of  the  optical  transfer  function  of  a  system 
for  different  angles  to  the  source.  The  two  parameters  r^  and 
9^  are  conjugate  pairs.  9^  looking  up  is  equivalent  to  r^ 
divided  by  the  path  length  looking  down  [Ref.  4]  and  vice 
versa.  A  more  formal  definition  is  that  9^  is  an  angular 
measure  of  spatial  coherence,  it  is  the  limiting  angle  for 
which  an  electromagnetic  wave  from  a  source  will  follow  the 
same  optical  path  length  to  a  receiver. 
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VACUUM 


Figure  1.  Comparison  of  LASER  Beam  Propagating 

Through  Vacuum  and  Atmosphere,  showing  r^ 


Optical  Paths  Optical  Path 


Figure  2.  Two  Sources  Showing  the  Effects  of 
the  Isoplanatic  Angle. 
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If  we  consider  two  paths  through  the  turbulence,  the 
isoplanatic  angle  relates  the  mutual  coherence  e'’  point  of  the 
B  field  between  the  two  paths.  Fried  [Ref.  9]  expresses  the 
isoplanatic  angle's  dependence  on  as, 

0^  -  2.905k“  C  “(z  z^/'^dz 

0  L  Jo  "  J  .  (6) 

where  z  is  the  altitude.  Note  the  spherical  weighting 

factor  that  emphasizes  turbulence  far  from  the  optical  system. 

The  preceding  paragraphs  clearly  show  the  importance  of 
the  parameter,  it  not  only  defines  the  measure  of 

turbulence,  it  also  determines  the  spatial  coherence  r^^  and 
the  isoplanatic  angle  6^.  The  problem  lies  in  that  high 
resolution  profiles  of  are  difficult  to  measure,  it 

requires  complex  detectors  and  optical  imaging  systems. 
Instead  we  can  define  a  temperature  structure  parameter 
similar  to  C^^  where, 

=  <(T2-Ti)  (7) 

which  can  be  measured  by  several  different  methods.  The 
fluctuations  of  the  index  of  refraction  are  due  to 
fluctuations  in  the  density  of  the  atmosphere  and  if  we  can 
assume  that  the  density  fluctuations  are  due  solely  to 
temperature  fluctuations,  then  is  related  to  by. 
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79x10' 

t2 


(8) 


where  P  is  the  atmospheric  pressure  in  millibars  and  T  is  the 
atmospheric  temperature  in  Kelvins  [Ref.  10].  The  assumption 
that  index  of  refraction  fluctuations  are  due  only  to 
temperature  fluctuations  and  that  humidity  fluctuations  are 
insignificant,  is  valid  when  the  Bowen  Ratio  B  (ratio  of 
sensible  heat  flux  to  latent  heat  flux)  is  greater  than  0.3 
[Refs.  11,12].  Below  this  value,  humidity  fluctuations  are 
significant. 


B .  SCALE  LENGTHS 

In  Kolomogorov's  definition  of  the  structure  function  he 
assumed  local  homogeneity  over  a  region  bounded  by  the  inner 
and  outer  scale  lengths  1^  and  L^.  These  scale  lengths  vary 
from  hundreds  of  meters  at  the  outer  scale  lengths  down  to 
millimeters  for  the  inner  scale  lengths  [Ref.  13] .  The  outer 
scale  length  is  the  size  of  the  turbulent  fluctuations  at  the 
onset,  while  at  the  inner  scale  viscosity  dissipates  the 
energy  of  turbulence  as  heat.  Kolomogorov  called  the  region 
between  1^  and  L,,  the  inertial  subrange  and  as  long  as  the 
distance  r,  in  Equations  (4)  and  (7),  is  within  this  inertial 
subrange  these  equations  are  valid.  Kolomogorov  expressed  the 
power  spectral  density  of  the  turbulence  in  this  region  by. 
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i 


(9) 


'!’(  K)  =  0.  033Cj^^K  ■ 

where  2ir  l;’  «  K  <<  2irlg‘\  Von  Karman  took  this  definition 
further  to  include  the  ranges  for  eddy  sizes  greater  than  L^, 

,-9/2  r  ,-"/« 

2^  .  3l,  ,  k2  I 

,  (10) 


r(ll/6)  i  1 

=  r(  1/3)  8 


1  + 


where  and  <S^>  is  the  variance  of  the  refractivity 

fluctuations  and  is  related  to  by, 


«  l.9<5„^K„^/^ 


(11) 


[Ref.  10].  Figure  3  shows  the  spectrum  of  turbulence  for  all 
scale  lengths,  the  inertial  subrange  shows  Kolomogorov's 
linear  description  of  turbulence  while  in  the  regions  above 
viscosity  effects  dominate  and  below  Von  Karman 's 

spectrum  defines  the  turbulence  [Ref.  14]. 

It  is  important  to  understand  the  effects  of  the 
turbulence  spectinm  in  all  three  regions  because  of  the  errors 
introduced  due  to  the  lack  of  correlation  in  the  regions  above 
and  below  the  inertial  subrange.  Care  must  be  taken  in 
choosing  the  correct  value  of  r  to  cover  the  scale  lengths 
inside  the  inertial  subrange.  Additionally  the  temporal 
frequency  of  the  turbulent  fluctuations  is  related  to  the  wave 
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number  K  and  the  wind  speed  moving  the  turbulence  past  the 
probes.  Therefore  to  measure  the  thin  transition  layers  in  the 
atmosphere  accurately  [Ref.  15]  the  errors  due  to  various 
scale  sizes  must  be  determined.  These  calculations  will  be 
carried  out  in  Chapter  4. 


SPECTRUM  OF  TURBULENCE 
FOR  nu.  St3<L£  LENGTHS 


Figure  3.  Spectrum  of  Turbulence  Over  All  Scale 

Lengths. 


C.  SYSTEM  REQDIRENENTS 

A  fast  response  time  temperature  probe  attached  to  a 
sounding  balloon  system  can  measure  the  vertical  profile  of 
the  temperature  structure  parameter.  The  design  of  the 
temperature  probe  system  provides  a  simple,  low  cost  method 


L2 


that  is  capable  of  resolving  the  thin  stratified  layers  of  a 
stable  atmosphere  and  determining  the  temperature  gradient, 
thickness,  and  turbulence  of  these  layers.  In  designing  the 
system  the  key  points  of  consideration  were  1)  what  type  of 
probe  geometry  and  thereby  what  type  of  circuitry  and  2)  what 
type  of  temperature  sensing  element  to  use. 

The  probe  system  could  use  a  differential  measurement  or 
a  single  point  probe.  The  advantage  of  the  single  point  probe 
is  that  it  can  measure  from  either  the  variance  of  the 
data,  knowing  the  balloon  ascent  rate,  or  by  analyzing  the 
power  spectral  density .However,  this  would  require  a  high  data 
rate  (several  samples  per  second) ,  which  the  radiosonde 
systems  used  do  not  have.  Although  a  differential  system  would 
not  have  as  high  a  vertical  resolution  it  has  the  advantage 
of  providing  partially  reduced  data  which  relaxes  the  need  for 
a  high  data  rate.  Therefore  a  differential  system  greatly 
simplifies  the  telemetry  needed  for  the  system  at  the  expense 
of  more  complexity  in  the  sensor  itself .The  resolution  for  the 
system  would  still  be  satisfactory,  about  2  meters  of  vertical 
resolution  for  a  balloon  with  an  ascent  rate  of  2  meters  per 
second. 

There  are  several  different  choices  for  the  probes.  They 
can  be  made  from  resistance  wires,  thermistors,  or 
thermocouples.  The  resistance  wire  is  simply  a  fine  wire,  such 
as  platinum  or  tungsten,  with  a  known  resistance  as  a  function 
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of  temperature.  The  problem  with  a  resistance  wire  for  this 
application  is  that  it  requires  some  means  of  self-calibration 
due  to  its  dependence  on  Ohm's  law  and  temperature.  A 
thermistor  is  a  small  semiconductor  device  which  changes  its 
resistance  as  a  function  of  temperature.  It  has  a  larger 
change  in  resistance  vs.  temperature  then  a  metal  wire, 
although  it  is  non-linear.  However,  thermistors  are  large 
compared  to  a  fine  wire,  having  a  larger  thermal  mass  than  a 
probe  made  from  a  fine  wire,  which  increases  the  response  time 
and  the  susceptibility  to  solar  heating  of  the  device. 
Additionally,  both  the  thermistor  and  the  resistance  wire 
r?.  quire  a  current  source  which  not  only  increases  the 
complexity  of  the  circuit  but  also  introduces  a  self  heating 
factor.  A  thermocouple  consists  of  two  fine  wires  made  of 
dissimilar  metals  welded  together,  which  produces  a  voltage 
difference  proportional  to  the  temperature.  For  an  in-depth 
discussion  of  thermocouples  and  the  thermoelectric  effect,  see 
Refs.  16  and  17.  Commercially  available  thermocouples  can  be 
made  of  very  fine  wire  (down  to  12.5  urn)  thereby  reducing  the 
thermal  mass  and  producing  a  faster  response  time.  A 
disadvantage  of  a  thermocouple  is  that  the  response  to 
temperature  is  small,  typically  40  ^V/°K.  This  places  severe 
requirements  on  low  noise  signal  processing.  An  operational 
requirement  for  this  type  of  device  would  be  the  knowledge  of 
the  mean  atmospheric  temperature  to  calculate  the  Seebeck 
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Coefficient,  which  is  the  derivative  of  the  thermal  emf  with 
respect  to  the  temperature.  Since  the  rawinsonde  system  used 
with  the  probe  provides  this  data,  it  is  easily  accomplished. 
Small  temperature  changes  of  1°  C  or  less,  are  expected  from 
the  two  probes  in  a  differential  system.  This  produces  a 
negligible  change  in  the  Seebeck  coefficient  therefore  no 
calibration  of  this  type  of  device  is  needed,  other  than 
knowing  the  gains  of  the  electronic  amplifiers. 
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III.  SYSTEM  DESIGN  AND  DEVELOPMENT 


A.  AMPLIFIER  CIRCUIT 

The  probe  system  consists  of  a  pair  of  thermocouples 
connected  in  series  and  held  rigidly  1  meter  apart  by  an 
aluminum  tube.  An  amplifier  circuit  requires  high  voltage  gain 
and  very  low  noise  to  be  able  to  discriminate  the  signal  from 
the  background  due  to  the  low  voltage  produced  by 
thermocouples  (on  the  order  of  /ivolts  per  degree)  .  It  also 
requires  an  analog  root  mean  square  device  to  calculate  <(T2- 

Ti)^. 

Since  the  purpose  or  the  system  is  to  measure  temperature 
differences  on  the  order  of  hundredths  and  thousandths  of  a 
degree  and  the  thermocouple  output  voltage  is  about  40 
microvolts  per  degree  C  the  circuit  must  have  a  very  high  gain 
with  ultra  low  self  noise.  Therefore  the  circuit  must  be 
carefully  designed  and  built  to  reduce  any  sources  of  noise 
wherever  possible,  by  methods  such  as  matching  resistor  values 
as  exact  as  possible  and  placing  them  as  close  as  possible  to 
each  other  to  reduce  the  thermal  drift.  Other  examples  include 
thermal  insulation  for  the  entire  circuit  and  RF  shielding  for 
the  circuit  as  well  as  the  probe.  The  most  critical  component 


of  the  circuit  is  the  ultra  low  noise,  precision,  high  speed 
op  amps  which  have  voltage  noises  less  than  those  of  50  ohm 
resistors. 

The  circuit  (designed  by  Prof.  Don  Walters  and  fabricated 
by  Dale  Galaro^'-'-cz)  is  a  low  noise,  wide  bandwidth 
Instrumentation  amplifier.  Three  operational  amplifiers 
produce  a  gain  of  10,000  and  a  combination  low  and  high  pass 
filter  with  a  gain  of  5  producing  a  total  voltage  gain  of 
50,000.  The  circuit  contains  two  Linear  Technologies  LT1028 
ultra  low  noise  precision  high  speed  op  amps  in  the 
instrumentation  amplifier  portion  of  the  circuit.  These  op 
amps  have  a  gain  bandwidth  product  of  75  MHz  and  a  self  noise 
of  0.8  5  nV/Hz^''^  at  the  .requencies  desired.  Figure  4  shows  the 
noise  and  frequency  characteristics  of  this  op  amp.  The  filter 
for  the  circuit  uses  an  LT1057  op  amp  whose  noise  and 
frequency  response  characteristics  are  shown  in  Figure  5. 
Based  on  the  frequency  response  curves  in  Figures  4  and  5 
[Ref.  IS]  a:.d  the  high  pass  filter  the  circuit  has  a  frequency 
response  from  0.16  Hz  to  200  Hz.  The  other  major  component  of 
the  circuit  is  an  Analog  Devices  AD637  high  precision,  wide 
band  RMS-DC  converter.  It  has  a  bandwidth  of  600KHz  at  lOOmV 
RMS  and  an  averaging  time  constant  of  25  msec/micro  F.  The 
entire  circuit  is  powered  by  two  9  volt  dry  cells.  Figure  6 
is  a  complete  schematic  diagram  of  the  circuit. 
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VOLIACE  NOISE  DENSITY  (nV/^Hz) 


Voltage  Noise  vs  Frequency 


0  1  10  100  1000 
FPEQUENCY  (Hz) 


Voltage  Gain  vs  Frequency 


c  01  0  1  1  10  If";  II  U">  IM  lOW  100M 

ff'l  C'JTNCY  I'*,-' 


Figure  4.  Noise  and  Gain  Characteristics  of  the 

LT1028  OP  AMP. 


Figure  5.  Noise  and  Gain  Characteristics  of  the 

LT1057  OP  AMP. 
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Figure  6.  Schematic  Diagram  of  Differentia 

Amplifier. 
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The  circuit  has  to  be  well  shielded  against  RF 
interference  due  to  the  high  gain  of  the  circuit  and  the  fact 
it  will  be  operated  in  the  atmosphere  where  it  is  highly 
susceptible  to  all  types  of  RF  signals.  Additionally  the 
thermocouple  wires  must  be  stretched  out  over  a  distance  of 
1  meter  and  will  act  as  an  antenna-  The  aluminum  tube  shields 
them  from  RFI  except  at  the  two  endpoints,  greatly  reducing 
extraneous  signals.  The  circuit  is  inclosed  in  a  styrofoam 
casing,  to  reduce  thermal  gradients  across  the  circuit,  which 
is  then  covered  in  aluminum  foil.  The  circuit  itself  has 
ferrite  beads  and  an  LC  filter,  at  the  inputs,  to  further 
shield  the  op  amps  from  RFI  and  the  entire  circuit  is  built 
on  a  ground  plane  which  has  been  grounded  with  the  foil 
shield. 

The  result  of  the  aluminum  foil  and  styrofoam  shielding 
and  use  of  ultra  low  noise  op  amps  is  an  amplifier  capable  of 
measuring  the  extremely  small  voltages  produced  by  the 
thermocouple  probes.  The  circuit  was  tested  in  the  anechoic 
chamber  in  the  basement  of  Spanagel  Hall  to  minimize  any 
temperature  fluctuations  and  then  measurements  of  the  circuits 
self  noise  were  taken.  Figure  7  shows  the  noise  spectrum 
measured  by  a  Hewlett  Packard  HP3561A  Dynamic  Signal  Analyzer. 
The  large  noise  spike  below  1  Hz  is  due  to  1/f  or  flicker 
noise  which  occurs  in  all  amplifiers,  due  in  large  part  to 
surface  leakage  of  transistors  [Ref.  19].  Although  the  noise 


normally  occurs  at  frequencies  up  to  lOOHz  it  has  been  reduced 
by  the  use  of  a  high  pass  filter.  The  noise  spike  at  6OH2  is 
due  to  AC  powered  equipment  operating  in  the  anechoic  chamber. 
Thx3  noise  should  vanish  when  the  system  is  used  in  the  field, 
since  it  is  powered  by  DC  batteries,  as  long  as  care  is  taken 
to  insure  it  is  not  near  a  large  AC  power  source.  The  plot 
clearly  shows  the  self  noise  output  of  the  circuit  is  well 
below  100  microVolts/Hz’''^  in  the  frequencies  of  interest. 
Since  the  circuit  gain  is  50,000  and  the  Seebeck  coefficient 
is  on  the  order  of  40  microvolts  per  degree  this  translates 
into  a  noise  induced  measurement  of  less  than  0 . 00005“c/Hz’^^. 


■  =  *.  :  -c  3rt;  954.83  S'OP:  ICC  -2 

Figure  7,  Noise  Spectrum  of  Amplifier  Circuit. 
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The  entire  system  was  set-up  and  run  in  the  anechoic 
chamber  with  caps  on  the  exposed  thermocouples  to  further 
reduce  the  signal.  After  running  for  about  2  hours  to  let  all 
air  currents  in  the  chamber  settle,  the  system  produced  the 
results  shown  in  Figure  8.  This  shows  the  noise  produced  by 
the  circuit  introduced  an  error  corresponding  to  a  of  less 
than  10'^,  which  is  two  orders  of  magnitude  lower  than  the 
lowest  needed  for  a  usable  probe  system.  These  results 
indicate  the  circuit  self  noise  is  well  below  that  which  would 
have  a  degrading  effect  on  the  results. 

The  entire  circuit  package  measures  3"  X  3"  X  3"  and 
weighs  less  than  six  ounces. 

B .  THERMOCOUPLE 

In  1821  Thomas  Seebeck  discovered  that  two  wires  of 
dissimilar  metals  joined  at  one  end  and  heated  produce  a 
voltage  difference  across  the  open  ends.  This  voltage  is  a 
function  of  the  junction  temperature  and  the  composition  of 
the  two  metals.  Since  that  time,  many  different  thermocouple 
types  have  been  produced  based  on  the  combination  of  the  two 
different  metals  used  and  having  different  thermoelectric  and 
physical  properties.  The  recjuirements  for  the  system  included 
use  in  the  atmosphere  from  the  surface  to  20  km,  therefore  the 
temperatures  vary  from  30°  C  to  -30°  C  with  as  large  a  Seebeck 
coefficient  as  possible.  The  Seebeck  coefficient  is  the  slope 
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Figure  8.  Noise  Output  of  Probe  System. 


Figure  9.  Photo  of  Circuit  and  Thermocouple  Probe 
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of  the  voltage  versus  temperature  curve  at  a  given 
temperature.  Based  on  these  reguirements  the  T  type  or  Copper- 
Constantan  thermocouple  was  selected.  An  alternative  choice 
would  be  the  E  type  or  Chromel-Constantan  thermocouple. 

The  Copper-Constantan  thermocouple  is  composed  of  a 
copper  wire  and  a  55%  copper  45%  nickel  wire.  It  has  a 
temperature  range  of  -200°  C  to  350°  C  and  is  suitable  for 
applications  where  moisture  is  present.  Table  1  shows  the 
thennoelectric  voltages  referenced  to  0°  C  for  a  Copper- 
Constantan  thermocouple  based  on  the  National  Bureau  of 
Standards  reference  tables. 

Table  1  shows  that  the  Seebeck  coefficient,  which  is  the 
unit  difference  in  voltage  for  each  temperature  change,  is  not 
linear  over  the  entire  temperature  range.  To  determine  the 
Seebeck  coefficient,  dV/dt,  valid  over  the  entire  temperature 
range  desired,  the  data  over  the  temperature  range  desired 
from  Table  1  was  plotted  and  then  a  polynomial  regression  was 
performed  to  find  the  equation  of  the  curve.  The  derivative 
of  the  curve  was  taken  to  find  the  Seebeck  coefficient.  A 
fifth  order  polynomial  fits  the  data  from  -100°  C  to  30°C.  The 
Seebeck  Coefficient  for  a  Copper-Constantan  thermocouple  is, 

dV/dt=3 . 87  07xl0'^+8 . 5348xl0'^t-3 . 3 13 5x10' V 

-2.774  32xl0’V-1.253xl0’”t'^,  (12) 
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where  t  is  in  °C  and  V  is  in  millivolts.  The  temperature 
difference  measured  by  the  probe,  is  the  output  voltage 
divided  by  the  circuit  gain  and  divided  by  dV/dt  determined 
at  the  specific  temperature.  For  a  1°  change  in  air 
temperature  the  average  change  in  dV/dt  is  less  than  0.2 
microvolts  per  degree  C  or  less  than  0.5%  change  in  dV/dt. 


1 


Figure  10.  Schematic  of  Thermocouple  Probe. 

The  probe  consists  of  two  0.00254  cm  diameter  Copper- 
Constantan  thermocouples  held  rigidly  1  meter  apart  by  an 
aluminum  tube.  This  wire  has  a  resistance  of  998.3  ohms  per 
double  meter  and  the  circuit  requires  a  low  input  impedance, 
therefore  it  is  necessary  to  cut  the  fine  thermocouple  wire 
short  and  solder  larger  wires  to  them  to  run  the  distance 
between  the  thermocouples  and  the  circuit  input.  Copper- 
Constantan  wires  0.0254  cm  in  diameter  were  used  which  have 
a  resistance  of  9.983  ohms  per  double  meter.  It  is  essential 
to  solder  the  copper  to  copper  and  constantan  to  constantan 
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to  insure  no  other  thermocouple  junctions  are  formed  which 
could  cancel  any  voltage  signal  generated.  The  resulting  probe 
yields  an  input  resistance  of  30  ohms. 


1/TflU  VS  SORT  of  WIND  SPEED 

FOR  Cu-Constantan  THERMOCOUPLE  .00254cm  DIRMETER 


(Hind  Speed  (m/sec)) 


C  .4  .8  1.2  1.6  2  2.4  2.8  3.2 

SORT  of  WIND  SPEe-D  Cm/sec  )^I/2 


Figure  11.  Response  Time  Versus  Wind  Speed  of  a 

.  00254  cm  Copper-Constantan  Thermocouple . 


The  reason  for  the  fine  wire  was  to  reduce  the  thermal 
mass  which  decreases  the  response  time  of  the  thermocouple  to 
temperature  change.  The  response  time  is  defined  as  the  time 
required  to  reach  63.2%  of  an  instantaneous  temperature 
change.  The  OMEGA  Handbook  [Ref.  20]  gives  the  response  time 
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for  a  0.00254  cm  thermocouple  as  0-05  sec  in  still  air  and 
0.004  sec  in  18.3  m/sec  air.  The  response  time  r  can  be 
related  to  the  velocity  by  the  equation, 

^  =  A  +  By/  V  ,  (13) 

T 

Where  A  and  B  are  constants  for  the  specific  thermocouple 
material  size  and  are  determined  by  the  response  times  given 
in  Ref.  20  (for  this  case  A=20  and  B=53.783  sec"’).  This 

results  in  a  response  time  verses  velocity  curve  shown  in 
Figure  11.  From  this  curve  the  response  time  t  §  5  m/sec  is 
.0071  sec.  The  interaction  of  the  frequency  response  and  the 
inner  scale  size  will  be  discussed  in  the  Results  section. 

C.  COMPUTER  AND  CODE 

The  data  acquisition  and  processing  portion  of  the  system 
consists  of  a  Hewlett  Packard  model  217  computer  with  a  20 
megabyte  hard  drive  and  2  megabytes  of  memory.  It  also 
contains  an  Infotek  BC203  Basic  compiler,  an  Infotek  AD200 
analog  to  digital  converter,  an  HP9133  floppy  disk  drive,  a 
monitor  and  printer.  The  computer  receives  data  from  the 
circuit  via  the  analog  to  digital  converter.  It  reduces  the 
data  and  displays  a  verses  time  (altitude)  profile. 

The  analog  to  digital  converter  is  the  primary  component 
for  the  data  acquisition,  however  it  introduces  errors.  The 


A  to  D  converter  has  a  small  DC  voltage  offset.  This  offset 
is  small  so  it  does  not  significantly  affect  the  results  until 
reaches  values  on  the  order  of  10’^,  as  is  discussed  in  the 
results  section.  Each  A  to  D  converter  has  a  different  offset 
so  each  time  a  different  computer  is  used  the  offset  must  be 
measured.  In  addition  the  noise  of  the  electronic  circuit 
produces  a  DC  offset.  To  find  the  total  DC  offset,  the  entire 
system  was  set  up  and  run  in  the  anechoic  chamber  to  insure 
no  external  signals  are  introduced.  The  software  has  a  feature 
for  inputting  the  offset  and  another  to  display  the  voltages 
directly  from  the  probe.  To  find  the  offset,  initially  set  the 
offset  to  zero  and  then  after  the  system  has  run  for  about  one 
hour,  to  settle  all  the  air  currents  generated  by  the  set-up, 
read  the  voltages  by  pressing  the  PRINT_RAW  key,  this  displays 
the  actual  voltage  offset  then  simply  reboot  the  system  and 
input  the  offset  read  from  the  raw  voltages.  The  DC  offset  for 
the  system  used  for  the  experiment  was  -5  millivolts. 

Appendix  A  summarizes  the  features  of  the  software  and 
lists  the  code  for  controlling  and  producing  output  from  the 
system.  It  is  almost  fully  automatic.  Once  it  has  started 
running,  the  only  corrections  to  be  made  are  to  update  the 
ambient  air  temperature  any  time  there  is  a  change  of  two  or 
three  degrees,  so  the  Seebeck  coefficient  can  be  recalculated. 
There  is  a  hard  function  key  for  updating  the  air  temperature. 
A  further  refinement  to  the  system  would  be  to  have  the 
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rawinsonde  temperature  measurement  device  automatically  update 
the  program  with  the  new  air  temperature  at  a  given  interval 
to  make  the  system  fully  automatic. 


D.  SOLAR  HEATING 

It  is  important  to  understand  the  effects  of  thermal 
radiation  on  the  probes  since  they  will  be  ascending  through 
the  atmosphere  both  during  the  day  and  at  night.  During  the 
day  radiative  heating  from  direct  and  reflected  solar 
radiation  will  heat  the  probes  above  ambient  air  temperature 
while  at  night  the  probes  will  be  cooled  due  to  Planck 
radiation  from  the  probes  to  space.  There  are  several  ways  in 
which  the  heating  or  cooling  of  the  probes  can  introduce 
errors  into  a  differential  measurement  device.  A  hot  wire 
anemometer  effect  can  introduce  errors.  Another  is  the 
difference  caused  by  one  probe  being  in  the  sun  and  the  other 
being  in  the  shade  or  at  night  where  one  probe  has  a  direct 
line  of  sight  to  the  earth  and  the  other  is  blocked  as  by  a 
cloud.  In  a  hot  wire  anemometer  effect  the  probes  are  heated 
above  the  ambient  temperature  and  then  velocity  fluctuations 
across  the  probes  vary  the  heat  transfer  rates  away  from  the 
probes,  creating  a  false  temperature  difference  [Ref.  21].  In 
the  sun-shade  effect,  as  the  probe  ascends  through  the 
atmosphere  one  end  may  come  into  shade  either  from  a  cloud  or 
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the  balloon  shadow.  The  probe  may  rotate  as  it  ascends  giving 
one  thermocouple  a  different  aspect  to  the  sun  ’-han  the  other, 
creating  the  same  sun/shade  effect  to  a  lesser  extent.  In 
either  case  one  thermocouple  will  receive  less  direct  solar 
radiation  than  the  other  thereby  changing  the  heat  flow  on  one 
thermocouple  but  not  the  other,  again  introducing  a  false 
difference.  The  same  applies  at  night  if  one  thermocouple  is 

shaded  from  the  earth  by  a  cloud  it  will  not  radiate  thermal 

energy  at  the  same  rate  as  the  other  will. 

The  net  heat  flow  to  or  from  a  body  in  the  atmosphere  is 
described  by  the  heat  transfer  equation, 

qnet=qs  +  qa  +  qt  +  '  (14) 

where  =  portion  of  direct  solar  radiation  absorbed 

q^  =  portion  of  the  atmospheric  radiation  absorbed 
q^  =  portion  of  terrestrial  radiation  absorbed 
q^  =  thermal  radiation  emitted  by  the  wire 
q,^  =  net  conduction  to  the  wire  from  the  atmosphere 

q^.  =  net  convection  to  the  wire  from  the  atmosphere 

[Ref.  22].  The  temperature  difference  can  be  determined  by 
first  setting  q^^^  equal  to  0.  The  equation  can  then  be  reduced 
into  the  heating  and  cooling  portions  by, 

^SD  ^SR  ^LE  ~  ^LW  '  (15) 
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where  Ejj,  =  heating  due  to  direct  solar  radiation 

Ejg  =  heating  due  to  solar  radiative  reflection  of 
the  atmosphere 

Ele  =  heating  due  to  long  wave  radiation  from  earth 
Em  =  cooling  due  to  long  wave  radiation  from  wire 
Ej.  =  convective  cooling 

[Ref.  23].  If  we  further  assume  the  thermocouple  to  be  a 
horizontally  oriented,  infinite  cylinder  with  the  top  half 
radiating  to  the  sky  and  the  bottom  half  radiating  to  the 
earth's  surface,  the  temperature  difference  between  the 
ambient  air  and  the  thermocouple  is  given  by, 


AT  = 


T  € 


+  R 


-  al' 


(16) 


where  h  =  average  convective  conductance 

€5  =  short  wave  emissivity  of  the  thermocouple 
=  long  wave  emissivity  of  the  thermocouple 
a  =  albedo 

R,  =  short  wave  incoming  radiation 
Rj  =  long  wave  radiation  from  the  earth's  surface 
R,  =  long  wave  atmospheric  radiation 
a  =  Stephan-Boltzmann  constant 
T  =  air  temperature 

where  h  is  defined  by  Kreith  [Ref.  22]  as, 
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where  D  =  wire  diameter 

V  =  wind  speed 

k  =  heat  conductivity  of  the  air 

V  =  kinematic  viscosity  of  the  air 

K&n  =  empirically  determined  dimensionless  constants 
based  on  the  Reynolds  number 

[Ref.  24]. 

The  Air  Force  Geophysics  Laboratory  [Ref.  23]  defines  the 
calculation  for  delta  T  in  a  similar  manner  based  on  the  same 
assumptions.  However,  when  comparing  the  two  forms  using  the 
same  parameters  (see  Figure  12)  there  is  a  considerable 
difference.  Review  of  both  treatments  shows  the  principle 
difference  lies  in  each  definition  of  the  convective 
conductance  h.  Campbell  [Ref.  24]  uses  Krieth's  form  of  b  (EQN 
17)  while  Brown  [Ref.  23]  assximes  an  average  value.  Another 
method  of  determining  h  is  with  the  Nusselt  number,  which  is 
a  dimensionless  number  used  in  describing  heat  transfer  and 
fluid  flows.  Kramers  [Ref.  25]  performed  extensive 
measurements  of  heat  transfer  on  spheres  and  cylinders,  from 
this  he  determined  the  Nusselt  number  to  be  a  function  of  the 
Reynolds  number  and  the  Prandtl  number,  another  dimensionless 
number  where. 
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Re  =  VDp  /  fi  , 


and  Pr  =  Cp^/k  , 


(18) 


where  D  =  diameter  of  the  cylinder 
V  =  velocity  of  the  fluid 
p  =  density  of  the  fluid 
p  =  dynamic  viscosity  of  the  fluid 
k  =  thermal  conductivity  of  the  fluid 
Cp=  specific  heat  of  the  fluid. 

From  all  the  available  data  he  showed  that  the  Nusselt  number 
for  a  cylinder  could  be  represented  by, 

Nu  =  0.91(Pr)°-^^  (Re)°-^®^  ,0.1<Re<50  , 

and  (19) 

Nu  =  0.60(Pr)°-^’ (Re)°-®°  ,50<Re<10,000. 

Based  on  the  Nusselt  number  the  convective  conductance  is, 

h  =  Nu*k/D.  (20) 

Campbell  included  some  experimental  results  in  his  paper. 
When  the  experimental  results  are  compared  with  Kreith's  and 
Kramers'  treatment  of  h  we  can  see  that  Kramers'  expression 
exactly  models  the  actual  data  (Figure  13)  .  Using  Equation  20 
for  h  in  Equation  16  corrects  the  differences  seen  in  Figure 
12,  therefore  the  lower  curve  in  Figure  12  is  the  correct 


model  for  solar  radiative  heating  in  the  atmosphere,  based  on 
the  results  obtained  using  the  experimental  data  of  Reference 
24 . 
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Figure  12.  Comparison  of  AFGL  Model  and  Campbell's 
Model  of  Solar  Heating  of  Thermocouple 
Wires. 


A  discrepancy  noted  in  Campbell's  calculations  was  the 
values  used  for  short (visible)  and  long(IR)  wave  emissivities. 
Table  2  gives  the  emissivities  used  and  the  actual 
emissivities  from  a  1986  edition  of  the  CRC  handbook.  The 
corrected  values  were  used  to  recalculate  the  solar  heating 
and  the  updated  results  are  shown  in  Figure  14 . 
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COMPRRISION  OF  CAMPBELL  TO  RFGL  h 


CAMPBELL (KRIETH)  _ AFGLCKRPMERS) 


EXPERIMENTAL 


VELOCITY  Cm/s] 


Figure  13.  Comparison  of  Convective  Conductance. 


Table  2 

EMISSIVITIES  OF  CU-CONSTANTAN 
AND  TUNGSTEN 


Visible 

IR 

1  Campbell's 

Cu-Constantan 

.25 

.5 

CRC  Handbook 

Cu-Constantan  .2  .03 


CRC  Handbook 

Tungsten  .5  .03 
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SOLRR  HERTING  OF  THERMOCOUPLE 


FOR  n  WIND  SPEED  OF  3  m/sac 


0  5  10  15  20  25  30 

ALTITUDE  CKm) 

Figure  14.  Corrected  Comparison  of  solar  heating. 


Based  on  Figure  14  it  is  now  possible  to  make  a 
determination  of  the  errors  introduced  by  solar  heating.  The 
structure  function  for  velocity  fluctuations  over  small  scale 
lengths  is, 

D(r)  =  3.83(er)^^^  ,  (21) 

where  e  is  the  dissipation  rate  [Ref.  26].  Actual  data  shown 
in  Reference  26  from  areas  of  highly  turbulent  velocity  fields 
indicates  the  average  dissipation  rate  is  on  the  order  of 
3x10'^  m^  sec'^.  For  the  probe  system  with  r  equal  to  1  meter 
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this  yields  velocity  fluctuations  on  the  order  of  0.06  m/sec. 
Figure  15  shows  the  solar  heating  errors  due  to  the  hot  wire 
anemometer  effect  with  velocity  fluctuations  of  this 
magnitude.  Figure  16  shows  the  differences  between  the  curves 
of  Figure  15.  It  indicates  the  temperature  differences  are 
negligible  (on  the  order  of  .001  degrees  or  less). 

The  shading  effect  can  be  determined  by  eliminating  the 
direct  solar  radiation  component  in  the  equation.  This  will 
show  the  maximum  error,  if  it  is  a  matter  of  the  probe 
changing  aspects  to  the  sun  the  errors  will  be  proportionally 
less.  If  direct  solar  radiation  is  completely  removed  there 
will  be  very  slight  heating  of  the  probe,  due  to  incoming 
terrestrial  radiation.  The  temperature  difference  between  the 
two  probes  will  be  approximately  equal  to  the  amount  of  solar 
heating  on  one  probe  as  seen  in  Figure  14.  This  indicates  a 
major  source  of  error  since  the  temperature  difference  at 
higher  altitudes  is  on  the  order  of  0.2  °C  and  the  balloon  can 
rotate  as  it  ascends. 

Based  on  these  calculations  if  some  method  is  devised  to 
limit  the  rotation  of  the  probe  to  eliminate  the  sun/shade 
effect,  the  probe  system  can  measure  values  accurately  up 
to  30  km  altitude  without  significant  errors. 

Reference  7  described  a  thermosonde  used  by  the  Air 
Force  Geophysics  Laboratory.  Data  measured  by  this  system 
indicates  an  order  of  magnitude  jump  in  the  values  of 
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Figure  15.  Hot  Wire  Anemometer  Effect  on  Copper 
Constantan  Thermocouples. 
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Figure  16.  Differences  Between  the  Curves  of  Figure 


just  after  sunrise  [Ref.  27].  This  increase  appears  to  be  an 
artifact  of  the  instrument  rather  than  actual  turbulent 
processes  due  to  the  fact  it  occurs  so  rapidly.  Since  it 
occurs  at  sunrise  a  logical  assumption  is  that  it  is  due  to 
solar  heating,  therefore  a  great  deal  of  time  has  been  spent 
in  determining  these  effects.  Figure  17  shows  the  hot  wire 
anemometer  effect  on  a  4  micron  tungsten  resistance  wire.  The 
values  for  the  emissivity  of  Tungsten  are  taken  from  Table  2. 
This  indicates  solar  heating  does  not  effect  the  measurements 
of  consequently  the  rise  in  the  value  of  at  sunrise  may 
be  actual . 

SOLRR  HEHTING  OF  4  MICRON  TUNGSTEN  WIRE 

_  FOf?  R  MEAN  WIND  SPEED  OP  3  m/sec 

*  ^  _  MEAN  VELOCITY  +/-  .06  m/sec  FLUCTUATION 
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Figure  17.  Hot  Wire  Anemometer  Effect  on  4  micron 

Tungsten  Wire. 
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IV.  RESULTS 


A.  EXPERIMENTAL  PROCEDURE 

The  experimental  measurements  served  two  purposes,  first 
they  were  carried  out  to  validate  the  probe  measurement  system 
and  second,  they  were  used  to  validate  the  measurement 
capabilities  of  the  acoustic  echosounder [Refs .  6,28].  The 

acoustic  echosounder  calculates  a  15  minute  time  averaged 
value  for  the  temperature  structure  parameter  as  a  function 
of  altitude,  however  uncertainties  in  the  antenna  beam  shape, 
side  lobes,  transducer  efficiencies  and  atmospheric 
attenuation  produce  uncertainties  in  the  absolute  value  of 
calculated  by  the  echosounder  [Ref.  28].  Therefore  independent 
verification  of  the  values  must  be  made  to  validate  the 
acoustic  echosounder.  If  both  values  agree  this  is  a  positive 
indication  that  both  systems  are  measuring  accurately. 

The  acoustic  echosounder  is  a  high  frequency  device  which 
analyzes  the  atmospheric  density  fluctuations  within  the  first 
200  meters  of  the  atmosphere.  The  echosounder  consists  of  a 
Hewlett  Packard  HP  217  computer  to  control  the  system  and 
acquire  and  reduce  the  data,  an  HP  3314A  function  generator 
which  produces  the  pulsed  signals,  an  amplifier  and  the  array 
of  speakers  which  acts  as  a  transmitter/receiver.  The  system 
operates  at  SKHz  and  produces  a  100  cycle  sinusoidal  burst  of 
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18  acoustic  watts.  The  antenna  array  consists  of  19 
piezoelectric  speakers  in  a  close-packed  hexagonal  array 
enclosed  in  a  55  gallon  plastic  trash  container  lined  with 
lead  and  foam  to  reduce  side  lobes  as  much  as  possible.  The 
minimum  range  of  the  device  is  approximately  six  meters  based 
on  the  recovery  or  "ring"  time  of  the  speakers  and  the  maximum 
range  is  about  200  meters  based  on  the  frequency  used.  Figure 
18  is  a  diagram  of  the  acoustic  echosounder  layout. 

A  comparison  test  was  run  with  the  probe  system  and  the 
echosounder  between  1300  and  2030  hours  local  time  on  3 
September  1988  on  the  roof  of  Spanagel  Hall  on  the  grounds  of 
the  Naval  Postgraduate  School  at  Monterey  CA.  The  acoustic 
echosounder  was  placed  on  the  sixth  level,  northwest  corner 
of  the  roof  while  the  probe  was  attached  to  a  rigid  pole  and 
extended  approximately  1.5  meters  off  the  seventh  level  of  the 
northwest  corner  of  the  roof.  In  this  position  the  probe  was 
approximately  9.2  meters  above  the  echosounder  array,  thereby 
being  outside  of  the  echosounder  blind  zone.  Due  to  the 
building  itself  and  heating  exhaust  vents  on  the  eastern  side 
of  the  building  it  was  necessary  to  monitor  the  wind  direction 
to  insure  the  prevailing  winds  were  not  passing  over  the 
building  and  picking  up  heat  from  the  exhausts,  which  would 
have  greatly  affected  the  data.  Therefore  wind  speed  and 
direction  as  well  as  temperature  and  humidity  measurements 
were  taken  every  15  minutes  to  update  temperature  and  humidity 
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information  for  the  two  systems  as  well  as  determining  if  the 
prevailing  winds  were  flowing  across  the  building  before 
passing  over  the  instruments,  corrupting  the  data.  During  the 
entire  experiment  the  wind  shifted  several  times  but  it  was 
always  from  Northwest  to  Southwest.  At  all  times  the  air  flow 
passed  over  the  instruments  before  passing  over  the  building. 
All  the  data  was  valid. 

Figures  19  through  21  represent  a  portion  of  the  lata 
collected  during  this  experiment.  In  each  of  the  figures  the 
upper  graph  is  the  echosounder  profile  of  the  atmosphere,  the 
central  graph  is  a  15  minute  time  averaged  profile  of  the 
atmosphere  based  on  the  data  collected  by  the  acoustic 
echosounder  and  the  bottom  graph  is  the  versus  time  plot 
measured  by  the  probe  system  at  an  altitude  of  9.2  meters 
above  the  echosounder.  The  dark  lines  below  about  6  meters  in 
the  upper  plot  and  the  discontinuities  below  6  meters  in  the 
center  plot  are  due  to  the  blind  zone  of  the  echosounder.  To 
compare  the  plots,  the  average  of  the  bottom  plot  was  compared 
with  the  value  at  an  altitude  of  9.2  meters  in  the  center 
pkot.  Figure  19  was  taken  early  in  the  afternoon  and  shows 
strong  convective  pluming,  which  causes  a  higher  temperature 
structure  parameter.  Figure  20  which  was  taken  closer  to  the 
neutral  event  shows  a  marked  decrease  in  the  turbulence  and 
a  corresponding  decrease  in  the  values  for  Figure  21 
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Figure  20.  Echosounder  Trace  and  Measurement  and 
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taken  during  the  neutral  event,  shows  virtually  no  turbulence 
and  a  much  lower  value  for  Appendix  D  contains  additional 
measurements  taken  during  this  experiment  to  show  the 
corresponding  increases  and  decreases  in  for  both  the 
echosounder  and  the  probe  system.  The  purpose  of  these 
measurements  was  not  to  actually  sample  the  atmospheric 
processes  at  this  location  but  to  make  a  quantitative 
comparison  test  between  the  acoustic  echosounder  and  the 
probe  system.  For  a  complete  description  of  the  atmospheric 
turbulence  measurements  and  processes  for  this  location  see 
Weingartner  [Ref.  6]. 

B.  SCALE  SIZE  ERRORS 

The  temperature  variations  in  a  turbulent  atmosphere  range 
in  size  from  millimeters  to  hundreds  of  meters.  Optical 
aberrations  are  primarily  caused  by  variations  the  size  of  a 
Fresnel  zone  therefore  with  laser  frequencies  and  path 
lengths  of  several  kilometers  the  important  scale  sizes  are 
on  the  order  of  several  centimeters. [Ref .  13]  With  a  frequency 
response  of  150  Hz  and  an  average  wind  speed  of  2-5  m/sec,  the 
system  is  limited  to  scale  sizes  greater  than  3  cm,  which  will 
introduce  a  small  amount  of  "inner  scale"  error,  since  minimum 
scale  sizes  are  on  the  order  of  millimeters,  but  if  used  for 
measurements  in  conjunction  with  laser  propagation  through  the 
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atmosphere  the  error  will  be  negligible.  Additionally,  the 
acoustic  echosounder  utilizes  the  smaller  scale  sizes  of 
approximately  3  cm,  thus  there  will  be  negligible  error 
introduced  by  this  in  a  comparison  test. 

To  find  the  outer  scale  length  errors  we  can  express  the 
structure  function  of  the  probe  system  with  limiting  scale 
lengths  by, 

Ut  “  s  i  (  k  r  /  2  ^ 


where  a  is  the  limiting  lower  frequency  (outer  scale) 

b  is  the  limiting  upper  frequency  (inner  scale) 

k  =  2n/A  where  a  is  the  actual  scale  length 

w  =  kV  where  V  is  the  wind  velocity 

r,  =  RC  ti-e  constant  of  the  high  pass  filter 
(the  upper  frequency  cutoff) 

=  frequency  response  of  the  probes 
(the  lower  frequency  cutoff) 

r  =  probe  separation  distance 

and  comparing  it  with  the  structure  function  over  all  scale 
lengths. 


D!0.x/=4j  sin 


r  .  2r^lk-5/3 
L 


(23) 


the  outer  scale  length  limiting  error  can  be  determined. 
Figure  22  is  a  graphical  representation  of  this  comparison. 
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showing  the  error  over  a  range  of  limiting  outer  scale 
lengths.  The  8%  limiting  error  on  the  low  end  of  Figure  22  is 
due  to  the  high  frequency  cutoff  of  the  circuit  and  as  the 
scale  size  decreases  the  larger  low  frequency  errors  of  the 
circuit  begin  to  dominate,  increasing  the  error.  The  limiting 
scale  length  for  the  experiment  can  be  determined  as  the 
height  above  ground,  which  was  approximately  30  meters.  From 
Figure  22  it  is  clear  the  error  introduced  due  to  finite  outer 
scale  lengths  is  approximately  12%.  The  design  of  the  acoustic 
echosounder  is  resistent  to  outer  scale  errors  therefore  the 
finite  inner  and  outer  scale  error  introduced  would  cause  the 
probe  system  to  record  measurements  approximately  12%  lower 
than  the  echosounder. 


f  7*1. z  I  ZC 


Figure  22.  %  Error  Induced  by  a  Limiting  Outer  Scale 

Length 
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The  acoustic  echosounder  is  susceptible  to  inner  scale 
errors  which  will  cause  it  to  read  higher  than  the  probe 
system.  The  inner  scale  length  is  inversely  proportional  to 
the  wind  speed  and  can  be  expresced  as, 

=  1  (24) 

Ochs  and  Hill  [Ref.  29]  made  extensive  measurements  of  the 
inner  scale  length,  based  on  their  results  and  the  mean  wind 
speed  of  6  m/sec  during  the  measurements,  the  approximate 
inner  scale  length  was  3  mm.  At  the  edge  of  the  inner  scale 
of  turbulence  there  is  a  bump  in  the  temperature  spectrum  due 
to  diffusion  as  it  enters  the  viscous-convective  range.  Figure 
23  illustrates  this  bump  showing  the  spatial  power  spectrum 
of  temperature  fluctuations  versus  the  scaled  wave  number 
KT) ,  which  is  the  wave  number  normalized  by  the  inner  scale 
length.  [Ref.  30]  Here  k  is  equal  to  27r/Scale  Length  and  r?  is 
equal  to  1^7.14  (for  air).  The  limiting  inner  scale  size  of 
the  echosounder  is  3.4  cm  (a/2  where  A  =  (340m/sec) (5kHz ) ) 
therefore  with  an  inner  scale  size  of  3mm  the  scaled  wave 
number  is  approximately  9x10'^.  Figure  23  shows  the  acoustic 
echosounder  will  read  approximately  5%  higher  than  the 
Kolomogorov  spectrum  and  therefore  5%  higher  than  the  probe 
system. 
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Figure  23.  Spatial  power  spectrum  of  temperature 
fluctuations  versus  scaled  wave  number  kt)  .  Solid  curve 
is  actual  model;  the  dashed  curve  is  Tatarskiis  model. 
[Ref.  30] 


C.  ANALYSIS  OF  DATA 

Figure  24  is  a  comparison  of  15  minute  time  averaged  data 
collected  from  the  echosounder  and  the  probe  system.  This  data 
was  taken  before  the  noise  measurements  and  discovery  of  the 
5  millivolt  DC  offset  error  in  the  A  to  D  converter  RMS  module 
combination.  Figure  25  shows  the  corrected  data  comparing  the 
two  systems.  The  data  clearly  shows  the  correlation  of  the  two 
systems  even  with  the  volatile  trends  of  the  turbulent 
fluctuations.  It  also  shows  a  decrease  in  the  temperature 
structure  parameter  leading  up  to  and  during  the  neutral 
event,  which  corresponds  with  the  actual  physical  processes 
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Figure  25.  Comparison  of  Data  From  Acoustic 

Echosounder  and  Temperature  Probe  After 
Correction  for  DC  Offset 
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TABLE  3 


c/  MEASUREMENTS 
(CORRECTED  FOR  DC  OFFSET) 


TIME r LOCALS 

ECHOSOUNDER 

PROBE 

%  DIFFERENCE 

1400 

. 0354184 

. 0368390 

3 

1415 

. 0316086 

.0391012 

23 

1430 

. 0358514 

.0335893 

6 

1445 

. 0500970 

. 0492531 

2 

1500 

. 0322137 

. 0321643 

0 

1515 

. 0490100 

. 0432243 

12 

1530 

. 0474621 

.0384319 

20 

1645 

. 0311984 

.0271304 

13 

1700 

. 0263566 

. 0350252 

32 

1715 

. 0361568 

. 0292298 

20 

1730 

. 0260007 

.0232678 

10 

1745 

.0217524 

.0274032 

25 

1800 

.0233976 

. 0200511 

14 

1815 

.0141386 

.0165680 

17 

1830 

. 0138539 

.0165680 

18 

1845 

.0132466 

.0184468 

39 

1900 

. 0129296 

.0145711 

15 

2000 

. 0054777 

.0077345 

40 

2015 

. 0058262 

. 0071837 

23 

AVERAGE  %  DIFFERENCE  17 

going  on  at  the  time.  With  the  approach  of  sunset,  at  1933 
local  time,  the  sun  heated  the  earth's  surface  to  a  lesser 
degree  thereby  reducing  the  temperature  difference  between  the 
earth's  surface  and  the  air,  which  in  turn  reduced  the 
temperature  structure  parameter  and  the  turbulence. 
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Table  3  contains  the  values  of  measured  by  each  of  the 
devices  and  corrected  for  the  offset,  it  indicates  an  average 
difference  of  17%,  with  the  probe  system  reading  lower.  The 
scale  length  errors  indicate  the  probe  should  read 
approximately  12%  lower  due  to  outer  scale  length  errors  and 
5%  lower  due  to  inner  scale  bump  errors.  Therefore  there  is 
no  significant  difference  between  the  probe  system  readings 
and  the  echosounder  calculations. 
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V.  CONCLUSIONS 


Independent  verification  of  values  measured  by  the 
acoustic  echosounder  is  important  [Ref.  28]  and  the 
differential  temperature  structure  parameter  probe  has 
provided  a  valuable  comparison  indicating  the  absolute  Cj^ 
values  of  both  the  echosounder  and  the  probe  are  valid.  Taking 
into  account  all  known  errors  there  is  no  significant 
difference  between  readings  of  the  echosounder  and  the  probe 
system,  which  is  an  extremely  good  indication  that  both 
systems  are  providing  valid  measurements.  Additionally  this 
thesis  demonstrated  that  solar  heating  of  the  probes  in  the 
atmosphere  does  not  appear  to  play  as  significant  a  role  as 
first  thought.  The  only  major  effect  solar  heating  has  on  the 
differential  system  is  when  one  probe  is  directly  illuminated 
by  solar  radiation  while  the  other  is  shaded. 

This  probe  system  has  many  applications  including  being 
placed  on  towers  to  calibrate  other  turbulence  measuring 
devices  as  well  as  being  attached  to  a  rawinsonde  system  and 
launched  to  measure  the  vertical  profile  of  When  used  in 
this  mode  it  can  measure  values  of  up  to  30  )cm  altitude 
accurately,  however  if  it  rotates  as  it  ascends  through  the 
atmosphere  the  sun/shade  effect  of  solar  heating  will 


adversely  affect  the  system.  Further  developments  to  the 
package,  such  as  addition  of  wind  vanes  on  the  probe  assembly 
which  will  not  affect  the  turbulent  flow  but  will  dampen  the 
rotation,  or  a  small  motor  with  a  flywheel  to  act  as  a 
gyroscopic  stabilizer  to  prevent  the  probe  from  rotating,  will 
eliminate  the  errors  induced  by  this  effect. 

Other  improvements  to  the  system  include  methods  to 
automatically  update  the  temperature  into  the  program  from 
the  balloon  systems  onboard  temperature  sensor.  Another 
improvement  would  be  to  increase  the  data  transmission  rate 
of  the  rawinsonde  telemetry  system  to  get  a  higher  resolution 
profile  of  the  thin  stratified  layers  of  the  turbulent 
atmosphere  or  even  possibly  having  the  system  transmit  an  AC 
signal  from  which  a  great  deal  more  information  can  be 
extracted  such  as  the  power  spectral  density  of  the 
turbulence . 
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APPENDIX  A 


PROBE  SYSTEM  SOFTWARE 

The  program  that  runs  the  system  is  called  "CTSQR".  It 
controls  the  probe  system,  collects  the  data,  reduces  it,  and 
then  displays  and  stores  it  for  further  analysis.  It  is  based 
on  the  same  program  that  controls  the  acoustic  echosounder. 
It  can  be  broken  down  into  several  sections.  The  first  section 
sets  up  the  system,  initializes  all  arrays,  creates  a  data 
file  which  can  store  up  to  eight  hours  of  data  and  sets  up  the 
function  keys  which  are  used  to  update  the  temperature  used 
in  calculating  the  Seebeck  Coefficient,  prints  out  raw  voltage 
data  or  ends  the  program  storing  what  has  been  collected.  The 
next  section  initializes  the  Infotek  AD200  analog-to-digital 
converter,  which  collects  the  data.  Now  that  the  system  is 
ready  to  collect  data  it  calculates  the  Seebeck  Coefficient 
based  on  the  information  input  at  start-up  or  updated  through 
the  function  key.  Next  it  collects  data  every  second  and 
averages  it  over  ten  seconds,  reduces  it  to  C^^  and  plots  it 
every  ten  seconds.  Every  15  minutes  it  prints  out  the  plot  and 
then  resumes  the  data  collection.  The  program  "CTREADER"  can 
take  the  data  file  generated  by  "CTSQR"  and  read  it  and 
calculate  15  minute  time  averaged  values  of  C^^. 
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integer  1 . J, Ji.K.Kstart .Kend.N.Nrec.Hr.O2<2000,4),Plotnun.Prlnt_key 
'  02  *  The  reduced  date  output  array  (2680,4) 

'  -  (Dev. Hr.  VOLTS.ftye,  CT  ) 

input  'ENTER  ftIR  TEHP  (DEGREES  C)“.T 


Se  t  cons  t  ant  8 
Pisc_addres8*** 

.700,0.6*  ' 

HPI8  addraea  of  disc 

Gflln-500e0 

1 

Anp 1 1 f 1 er  gain 

Maxrec-rPPB 

( 

*  records  in  output  file 

NdIcI-900 

( 

f  point s  plotted 

Piet  nun-0 

Pn  1  nt  _1<  ey  •  1 

1 

print  raw  data  if  ^0 

R-1 .0 

1 

Probe  separation  (n) 

Scele-10000 

‘  Olac 

storage  scale  factor 

’  The  Eouation  Per  The  Seebect^  Coefficient 

See-3 .e7?7F-2*e.5340E'S*T-3.313SE-7*T-2 
Pec) --7 .77432E-9*T-3-1 .253E'1 1*T'4 
Seebeefc  •<  See»Eec)'  'M.E-3 
PRINT  Secbeck 

I 

P__cne^(hird*R'(  1./3.  ) 

t 

'SE^  TiMECftTE 

input  *00  YOU  WANT  TO  RESET  THE  CLOCK  (Y  OR  N)7',0t 
IF  0***Y*  THEN 

INPUT  ‘ETITEP  **00  HHtt  YYYY"  (Local  Tine)*. Dale* 
input  'ENTER  "HR  MIN  SC**  (Local  Tine)*. Tine* 

SET  riHEOftrS  DftTE(Oete«)*TlM£(Tine*) 

POINT  oftie*( timedate ). tim£»( timedate > 

Titert-TIHCPATE 
T?-Tst-^rl  MOO  8G400 
FNP  IF 

INPUT  'INPUT  site  flAME'.Slte* 

INPUT  TNfFP  the  A-0  converter  OFFSET'-. ??5  FOR  HP  2l7)'.?n'*o 
IflPur  'FNTFP  THE  LOUEST  OECftCE  FOR  THE  PLOT  INORMftL  USE  -5)'.Ynin 

'  S*l  un  (Ne  d*te  reduction  outnut  file 

lNr"T  'FNTFP  nEOOCED  DATA  OUTPUT  FRE  NftME’.ElleJl 
Pilel5-Pliel*iDl  s*‘_cddre5B» 

INP'.'T  MSr  ENTPy  IN  reduced  OUTPUT  FILE7  (yES  OR  N0)*.Q» 

IF  Of-’NO*  THEN  GOTO  Oldfile 

M-wflle  CPEftTF  POAT  F  1 1  e  I » .  1  . 23040  1  2  BYTES  FRE  SIZE  6  HOURS  OF  DATA 

ASSIGN  AF  I  lei  TO  F  I  Ult 
N.-CC-?  'I  OF  entries  in  the  OUTPUT  FR£ 

GOTO  SeluD 

Oldfile  ASSIGN  pFilel  TOFilel* 

ENTER  ?Fi lel .D2( •  ) 

ASSIGN  fiFllel  10  Fliel*  tGO  TO  START  OF  FRE 

Nrec-P2<J.l)  '  THE  OLD  t  OF  ENTRIES  IN  THE  OUTPUT  FRE 
'etip  '  Set  uD  the  deta  reduction  and  olol  fornet 
erPUT  KPO. 'SCRATCH  KEVE'i  «  Clear  keys 

C^NIROL  2.7,1  I  Select  user  nenu  I 

ON  ypi  1  lAEEL  ‘PRINT  RAW*  GOTO  Prinl^rew 
OH  KFf  0  IPOFL  'Quit*  GOTO  Quit 
ON  KEY  2  LAPEL  •UPDATE  TEMP*  GOTO  UnHate.tenp 

NnrlPt-TIMEonrE  MOO  66400  MOO  3600  DI'J  15 
CAI L  P I c t • e (up ' Np 1 ct , Si te* . Ynin ) 

pall  Init  «d2?0  (  Iniljelira  A-D 

'  OUTPUT  KEI’  'L*.  ’  Turn  on  grephtce 

'  Begin  the  nein  data  acqulaillon  loop 

I F  Nr-C  'r  er 
:  e  ,  1  ■'  ' 
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76? 

770 

790 


Vol i »Q-0 
Store_dele-0 
FOR  1-1  TO  10 


790  Sync'  '  5ynchronl2B  data  collection 
800  T1-1NT(TIMEDATE  HOO  86400) 

eie  IF  T1<T0  THEN  T0-T0-eG<00 

920  IF  (TI~T0)<1  THEN  60T0  Sync 

930  T0-T1 

840  Read.ad'  • 


with  the  By8le«  clock 


6^0  CALL  Adln(UoltaQe,7ero)  t  Read  Inttiak  A-0 

860  IF  Print_key>0  THEN 

070  PRINT  USING  M0D.OOOO-:UoltaQe 

880  GOTO  Readied 

890  END  IF 

900  OeltsQ-Uoltaq+Uoltaoe^Uoltage  I  Average  voHage*2 

910  Npolnt-lNT( Ti  HOO  3G00  MOO  Nplot) 

920  IF  Moolnt <NDOtnl_old  THEN  Store.data-l 

930  IF  Npolnt_old>0  AND  Store_dBta«0  THEN 

940  I  Plot  the  data 


350  C t aor«(Uo It  age/ ( Gal n«Seebeck«R_one_ third ) )*2- 

960  Lgclsor-LGItCtaqr ) 

970  I  MOVE  Npolnl^old.Vol ta_old 

9P0  '  DRAW  Npclnt  .UoMage  i  Plot  Uoltage 

990  MCUE  Npo  1  ni _ol d, Lgc t BQr__ol d 

10®0  DRAW  Npoint .Lgctsqr 

1010  ENP  IF 

1020  Vc 1 t 5_o 1 d-Uol t age 

1030  Loct5Qr_old-Loct5Qr 

1040  Mpc 1 nt_o 1 d-Npol nt 

1050  IF  5tore_data-l  THEN  Npolnt^old-G 

1060  NE^T  I 

1070  Vcll«_avo«$ClR(yoltsQ/10) 

1080  Ct-VoUs_avg/(6ain«Seeback*R  cne  third) 

1090  DISP  •  • 

1100  n-TlMEPATE 

1 110  Payt-PATE*?  T1  ) 

I  120  Tinet-T1ME*( T1  » 

1130  • 

1M0  '  CT  SCUAPEO  OATA  REDUCTION  SECTION  ••• 

I  1S0  ' 

1 160  vrf-Oe, 1i 8.  I  1  ] 

1170  Cay-( ( PATE<0ev*  )-0ATEt • 1  JAN  ’GYr*))  OIV  0B40e)*J 
1180' 

1  1  90  Mre-*^lrec*  1 


1200  '  2At.C''l.ATE  OECIHAL  HOURS 

1210 

1220  M-uri-'.'AL'  Ttn  .?  I  'MUAL(  T$(  4.5)  )*VAl  t  T»I  7 . 8  )  )/60  )  '60 

1220  Hf -1000*Hcurs  '  Mete  that  the  HP  rounds 

12^^  rniMT  *Oeccrd  f.Mrec;*  Collected  •.Oayt;’  ^Tjnef 

i:r(»  coiMi  -  • 

1  2F0  ALPHA  prF 

ir**0  '  5ET  I'P  OUTPUT  APr.^v 

1280  fl!-Nre-*l  «FPR  PPTIC-N  °ASE  I 

12^0  PT' 1 . 1  )-M-ec 

13^0  P2'»il  .1  '-Pay 

12  10  02'  rii  '-Mr 

12.0  Dr'fJl,3)*‘.'olt*_avg»Scele 


l2-»0  P?'M1  .4  )-CfScale  i  IntenMtv 

1.’40  PRINT  N1  .C2(Nt  .  1  '.0?<N1  .r  ):P2TMl  ,2  );D7(M1  .  4  ) 

1260  Slc^e.data  »  U*- )  t  a  output  every  NPLOT  eecends 

I.TF0  IF  Stcre  dala-j  THEN 

1 27P  ^ \ cr e_da  t  e-0 

I  ■'P0  PRINT 

12R0  PICO  'WRITING  PFDUCEO  OUTPUT” 

1400  Output  fr j l e l .02( • > 

14  10  a<5c  ip.,  ftp  j  Igj  TO  Fi  lelt 

1470  Plrlntin-PlctnuM*! 

H’p  IF  Plctnun  MCD  2-1  THEM 

1 140  PRINTER  IS  701 

PPirtT  • 
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!4G?  PRINTER  IS  i 

EMO  IP 

HP?  D’JMP  GRAPHICS  I7®1  1  Ou«p  screen  to  printer 

149?  Npclnt-? 

IGP?  DISP  ■  ■ 

ISie  CALL  PlotseiuptNplot.Sltfif  ) 

152?  EMO  IF 

IS??  End_«Hi le  END  UHILE 

154?  1 

155?  Updole_tenp-  •  Updates  SeebecP  Coefficient  With  New  Air  Tenp 

ISE?  INPUT  -ENTER  NEU  AIR  TEMPERATUREt DEGREES  O'.T 

157?  See-3.87?7E-2+8.5340E-5*T'3.3135E-7*T-2 

J5P?  Peck--2. 7743?E-9*T*3- I .7S3E-1 J»T'4 

159?  Seebeck“<See*0eck)*l.£-3 

16??  GOTO  Start_lo 

IBl?  Prlnl_ray  >  TOGGLE  THE  PRINT  FLA6 

1B2?  Print_k  ev"-Pr lnt_key 

IE??  GOTO  5tart_lo 

164?  Quit  *  Urite  reduced  data  output  file 
165?  FOR  t-1  TO  N1 

165?  PRINT  I :02t 1 . 1  ).02t I .2  ).02< 1 ,3);D2t I .4  1;D2( 1 .5  );02{ 1 .6  ) 

1E7?  HF^T  I 

166?  OUTPUT  PFt lei ;D2< • » 

169?  PRINT  'DATA  FILE  HAS  BEEN  STORED  UNDER  NAHE'.Filel* 

17??  PEEP 

171?  PEEP 

172?  A55IC-M  CFl  le!  TO  • 

173?  CTCP 
174?  END 

175?  SUP  P 1 c t set up< Np lot , S 1 t e$ , Y«ln ) 

I7E?  Ynex-1 

177?  C-INIT 

17P?  6PAPMICS  ON 

179?  LIME  type  1 

le??  UIFUFCRT  15.120,15.0? 

1610  WINDOW  e.Nplot .Ynln.Ynax 

162?  AveS  60. .S.?.VBin.5.2 

16!?  CLIP  OFF 

1P4?  CSI2E  4.  .6 

165?  LCPC-  6 

166?  '  P-au  Leg  Y  A<ie 

167?  *  for  recade*Y«in  TO  Yrta'< 

1660  '  FOR  l'nits-1  TO  1 -tSM  Decede<Ywe'<  » 

169?  '  Y«Dec ede *L6T ( L'nl  1 3  ) 

la??  '  MDUE  ?.Y 

191?  '  DRAW  Nplct.Y 

192?  '  MEYT  Units 

192?  '  NEYT  Decade 

194?  I  Label  hrrijcntal  axis 

195?  TI.TIMEDATE  HDD  66400 

19S?  Hrs-n  Diy  3600 

197?  tr-Tl  MOO  360? 

1°6?  Mm-T?  OIU  G? 

199?  Otrhr-Mio  DIU  15 

2???  FOR  M-0  TO  NpJot  STEP  300 

771?  MOVE  M.Y«ln-.?6 

202?  Olrn{n-Otrhr» IS*<M/300 )*5 

2?’?  IF  OtrRln-60  THEN 

204?  Otrnln-0 

205?  Mrs-Mre^l 

206?  FND  IF 

??7?  l.APFL  USING  ‘OD,A.2Z’.Hrs;’  ■;0trnln 

2?P?  NE'^I  M 

209?  Move  NDlot/2.Ynln-.0 

71??  LAPEL  -Tine  tLocal  )’ 

211?  *  Label  Ordinate 

712?  LCPG  6 

217?  FOP  M-Ynln  TO  Yns- 

214?  LOP.r*  B 

21*^7  CSIZE  4 

2IF?  MOVE  “Nplct/25.3.M 

717?  LAPEL  USING  -I.K'.-l?' 

216?  CSIZE  2 

2  19?  l.OPG  I 

2  7??  MP'.T  -Nolct'76.M 

Zri?  LAPEL  USING  ’I.K*  M 

:::?  ne'^t  m 
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Mt-'f  M 

hOVE  Np  lol /2  .  Ynln- .  0 
LAPEL  ‘Tine  tLocal)' 

Zn?  >  Label  Crdlnala 
212?  L0R6  e 

21??  FOR  M-YJiiln  TO  Ynas 
2M?  LORE  0 

215?  CSIZE  * 

2160  MC'.'E  -NdIoI/23.3.H 

2170  LAPEL  USING 

ne?  C5I2E  2 

210?  L0R6  1 

22??  HO'.'E  “Mplet/Ze.H 

221?  LAPEL  USING  't.K" :h 

222?  MEYT  M 

22??  LOIR  PI/2 

22<?  LORG  6 

22=?  CSIZE  4 

226?  MOUE  'Np I c* /7, <  )/2 

227?  LAPEL  -LI’S  OF  CT‘2" 

226?  ITil  le  the  plot 
229?  LOIR  ? 

23??  LORE  4 

2?  I?  f’CUE  MclcW2.Y«a«*I 

232?  LAPEL  PATEKTlMEDATEli'  ‘;SMet 

23'’0  Cl  IP  ON 

234?  SUPEND 

235?  sue  Adln'VcMage.Zerc) 

236?  '  26  APR  ISPB-  DLU 

237?  '  IMF0TE><  A-0  input  routine  set  uo  for  internal  trlOOer 

23??  ’  and  average  40  points  over  three  60  Hi  cycles 

239?  ' 

2  4??  IMTEC-ER  I  .Nooints,Ad_daia<  1-40) 

241?  OIM  Select»[2?I 

2<r?  Ad_5el _ccde»] 7 

2  4  3?  1 nt  eni 1  I  V*? 

244?  Mpe i n i s ■ 4? 

245?  Ccuntt-UALtlHpCints ) 

246?  Scele-5.?/lNpoints»2?4?. ) 

247?  Sldrv-0 

24??  Del te_t ine*** 1250???*  1  AO  interval  between  senples  in  nsec 

2450  SeUctf-'eelect  Isl  end* 

25??  GCSue  Read_od 

2?l?  FOR  1*1  70  Npoints 

7^2?  '/c  1 1  age  •Vo  I  tegc*Ad_delo<  I  ) 

753?  fiEYT  I 

75  4?  Vc 1 1  oge«<  Vo  1 1  age-Zero  )*Sce I e 

255?  GOTO  Subend 

756?  Reed.ed  '  Reed  the  Infotek  A-0 

257?  i  Initleliae  the  A-0 

25P?  OUTRL'T  Ad_se  l^code  • 'RESET" ,  •  I nterna  1  *  , ‘count  * ,  Coun I  f 

25=?  OUTPUT  Ad_eel_code ; M  me' .Del  tontine*. 'deleyon' ,5elec t« 

26??  OUTPUT  Ad. sol.code: -STATUS’ 

2610  ENTER  Ad.sel^code iReen* 

2620  IF  Peso*-’ - them 

2630  ENTER  Ad.sel^codc  USING  M.W  ;Ad.data(  •  > 

2640  OUTPUT  Ad_seT_code; -STATUS' 

265?  ENTER  Ad.ee l.code : Reso» 

266?  IF  Resp»<>- - '  THEM 

267?  PRINT  "ERROR  DURING  SAMPLING  •  •;Re«p* 

269?  END  IF 

269?  ELSE 

2'’??  PPINT  ’ERROR  DURING  A-0  INITIALIZATION  -  "vRaept 

7710  ENO  IF 

27??  RETURN 
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2730  Subend'  * 

2710  SUBEHD 

:7r.O  SUB  Init_ad2P» 

2760  I  initialize  0D_200 
2770  Code-17 

2780  Ounny-BEAOIOl  Code , 3 ) 

2790  URITEIO  Code, 010 

2800  CONTROL  Code, 01 1 

2810  5UBENO 
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APPENDIX  B 


SOLAR  HEATING  PROGRAM 


13 
'3 
33 
13 
S3 
S3 
73 
89 
99 
1 39 
[  '  9 
129 
1 23 
119 
1E3 
1S9 
179 
189 
199 
290 
219 
229 
229 
210 
2S9 
2S9 
2'’9 
299 
290 

tttix 

219 
229 
229 
319 
259 
269 
379 
399 
299 
199 
119 
129 
129 
119 
1S9 
169 
♦  79 
199 
199 
.  599 
519 
5  1 1 
512 
512 
511 
529 
529 
519 
cr  •> 


9E-ST0RE  VurPLSTl  .730,1, 9‘ 

,  nr  T  )  ooo  “Rn 

THIS  P9CSP1W  °1-3TS  the  CHAM6E  IM  THE  SE'-T-l  T  '  TEHPEPRT'JPE  OISSEPEMCE 
ggTijcr»j  TH^  THE9.HCC2'JPL£  RMC  ^HE  .IIP  TE“8  CUE  TO  cm  HE^T^*J5> 
ijjTH  IMCPEI^SE  111  RL^IT'JEE  9.1155C  3M  SCCC’S  C*LC*JI_.1T  lOMS  IM  l^FSL  P'J9 
O-TEO  22  S5S  1991  .‘IfiD  CPHPSELL’S  UGPH  FRGH  GOT  ;9S9 


1  THE  FGLLGUIHS  IS  THE  LIST  GF  '.'P..91P2L5S 

i 

J'4.  19 
£*.  M 
E53t'.:5 

C^-«0  T 

Fr-2. 

Ec^iw'.S 
£=3l'.£ 

Ro-.e:: 
fl^-.399 
Rg-.aiS 
Albff5c'.:£ 


‘  *^«han:cai  scuivalent  zf  hftjt  Cjl"-1  Sec' 

‘  Soier  Careter.t  ^W  c?t‘-2) 

*  VieiSle  yeveie-Sth  A8£0RPTI«*Ty  zf  t-e  TM£®“CCC'JPL£ 
'  Shor^  Ueve  Enieeivlty  fzr  CA?1PSELL’£  Celcjle*  icr.a 
'  Pcrf%  Pester  fzr  Cir^ct  fte^ietisn 
'  Per."*  Fester  ?cr  Reflected  Redtetisn 
'  Long  yeve  EffISSiVlTY  cf  t.*^e  TH£R?fCOC'JPL£ 

>  Lsrg  yeve  Eflteeiv;*/  fer  CAftPBELL’S  Ceisulet  tsne 


'  Shsrt  Weve  Ircsr^irg  fledietlcn  tCel  cn"-2  3ec*-l) 

'  Lsng  Weve  At5»oecheric  Redietlon  (Cei  cr^‘*2  »es*-l  > 

‘  Lsrg  yeve  Redletlon  frzn  3rsun<l  'r>!  sc»“-2  3ec‘-l) 

'  Reflectlsr*  of  the  Eerth 
'  Stefen-BoHci^enn  Cenetent  tW  cw*-2  K‘-*> 

Sigcel  . 25E99E-12  »  Stefen^ol  tc»enf>  Conatent  (Cel  cn*-2  K*-4> 
O«*.00992S4  I  Prebe  Oleweter  in  neter^ 

Oci^'*. 00254  ‘  Probe  Oie^^ter  in  senti»^etero 

Csr^'i'ert '•i  .  1 50025^2  ‘  Corweroion  Festor  fzr  Ther«el  Cc^dus  *  1  •  J  t / 

Retell . -iSSE-S  ‘  Conatent  for  Oeternmirg  hu  (kg‘-l  r.‘-i 
Suth-«110.4  *  SutherleriJ'a  Csriatent  for  ffe  (X  ; 

Pr-.?i4  ‘  OlJ^naionleee  Prendtl  Nvf'-.oer  fc**  4:r 

iRe  ••  C:»eraisnieoe  Reynolds  Nunoer  fzr  air 

'fkid  ^  Oinenaionieoe  Nueaeit  Nunber  for  A;r 

IT  ••  Air  Tenoereture  'Kelvira? 

iRcw  ^  Air  Oenaity  eo  Function  of  Altitude  ?kg/n*2> 

(rtu  •  Oynemc  Vlacoe:ty  of  Atr  {N-aec/n*2^ 

‘Hu  ^  Kirenetic  Viacoeity  of  Air  (n*2/aec> 

IK  ••  Thernel  Conduct  :*  1 1>*  (Cel  >ec‘-l  cn‘-2  (C/cn.'“-l> 

‘H  Convectivr?  Meet  Trenofer  Coefficient  (f^ALT;) 

I  (Cei  cn‘-2  aec*-l  K'-l > 

»F  •  Percent  go  of  Soler  Redietion  Reechirg  5iven  Altitude 

'So  '*  Meet  ®'lux  fron  Eerth  et  5lven  Altitude  (W  cn"-2.' 

(Ait  •  The  ftiven  Altitude  (Kn> 

"j  -  yind  Soeod  (n/aec.' 

Eneiricei  conotent  beoed  cn  ReynoJda  Hunter  fren  Krieth 
E^»oiricei  cenotent  beoed  on  Reyncida  Hunter  fren  Krieth 
Tenoereture  Olf  Terence  between  TC  j,r:d  Air  for  SOCO 
Tenoereturo  Difference  between  TC  end  Air  for  CAAPSELL 
9s th  in  Oegreeo  C 


iKr  - 
(H 

’Deit' 

'Delti- 

i 
\ 

INPUT  'WHAT  WIHO  SPEED  CO  YOU  WAMT*,V 


I  pi^nr  CC'^’jP 

INPUT  ’INPl'T  1  erp  PLOTTER  CP  2  PCR  CRT’.C 
IF  Q'l  Them  olC^TER  IS  707,-mp5L* 

1  BE  SURE  RLCTTER  CIP  SWITCHES  ARE  P°CP£RLY  SET  le.  SWI'^CH  1.2,2  IN  PC5IT  1 
IF  2'2  I'MEn  PLCTTER  IS  CRT INTERNAL  ’ 

I  SIHI’’ 

I  cpdou^-^c  '2^ 
t  » --C  •  ^OC  ' 


.E.  ;:a. 


'.90 
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5E2  ‘J!MCOU  1 

c“2  AXES  ^  2^  **  3  5  2 

522  CLIP  CP^ 

c ort  C5 I Z£  4  .  S 

ex*  <  '^QC  g 

si2  ‘  label  horizontal  axis 

£:0  FCR  H-e  "^0  30  STEP  5 

£23  ?1CV£  rt,-.04 

EAC  LABEL  H 

£5?  •‘lEXT  H 

SEC  HOVE  15,- . 1 

£“'3  LABEL  “AL’’IT‘JCE 

ES3  '  LABEL  VERTICAL  AXIS 

can  LOPS  S 

"^32  ROR  TO  1  STE°  .  1 
’’13  “OVE  -  5  H 

n-  ft  LA®E'  H 

M^XT  “ 

740  LOIR  PI/2 

'’53  L0R5  S 

-fcn  ksCitic  _  c  c 

label  *05'  T  'Tx* 

•^00  I  TTTIC-  pir\T 

^  om  I  r>  r  Q  •” 

.  ^  W  W  k  ■  '  w 

93»  LGRS  * 

o '  0  ^  *  ’ 

203  LABEL  "SCLAR  HEATIN&  OF  THERMCCC'JRLE * 

932  f^fC'.'E  IE,  I 
510  CSIZE  3 

353  label  ‘FOR  A  'JINC  SPEED  OF’ j  V; 

SS3  CSIZE  * 

073  I 

352  ^  C^lcul^tisns 

333  ' 

323  AMci2'^ 

313  Ce-lSsld*^ 

923  C^iilold-*2 

523  FCR  Alt-a  TO  29 

349  ^  I  w  Jssunir.g  the-  relatlcnshic  for  F  and  Sa  ara  l;r.*jr  jrt  Alt;*.ud 

952  ‘  Tha  Valjaa  wara  takan  frew  BROWN  and  6CC0 

552  F'.SXAlt  v91.2333  > 

9'’3  Sa*.  9E£-<  Al*.  •  .90141925  ) 
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